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Abstract 
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Introduction 



The recent progress in calculating correlation functions of quantum solvable 
models is based on the fact that they are governed by classical integrable dif- 
ferential equations. That the language of classical differential equations is quite 
natural for the description of quantum correlation functions was realized a time 
ago |], |], f§] . The idea of the approach suggested in |5], H, [7| is to consider the 
Fredholm determinant in the representation for a correlation function of a quan- 
tum integrable model as a tau-function for a classical integrable system (see also 
the book || where the results for the simplest model of one- dimensional impene- 
trable bosons are reviewed). The necessary first step (which is also of interest by 
itself) in this approach is to represent the correlation function as the Fredholm 
determinant of a linear integral operator of a special kind (the "integrable" in- 
tegral operator). It is to be mentioned that the first determinant representation 
of this kind was given in || [H| for the equal-time temperature correlators of the 
impenetrable bosons in one space dimension. 

We consider the exactly solvable one-dimensional model of the nonrelativistic 
two-component Bose and Fermi gases with the infinitely large coupling constant, 
c = oo. This model for any values of c was solved in the paper of Yang [jllj where 
the eigenvectors of the Hamiltonian were constructed using relations known now 
as the famous Yang-Baxter relations. The general approach to the solution of the 
model with n internal degrees of freedom was given in paper see also [13 



(for the solution in the frame of the algebraic Bethe Ansatz see paper [Q). The 
Fredholm determinant representation for the equal-time temperature correlators 
of the Fermi gas with c = oo was derived in paper jnjj where also the long distance 
asymptotics was constructed. 

In our paper the Fredholm determinant representations for the time dependent 
temperature correlation functions are obtained, for both Bose and Fermi gases 
with c = oo. To do this, we generalize the method used earlier for the one- 
component models [jTBL [T7|, [TS|, |T9fl . A presentation of our results was given in 
paper JSJ. 

The contents of this paper is as follows. In Section 1, the explicit expressions 
for eigenstates (the nested Bethe eigenvectors) are given. At c = oo the auxiliary 
lattice problem for the nested Bethe Ansatz is reduced to the problem of finding 
eigenvectors for the cyclic shift operator and one can use the eigenvectors of the 
Heisenberg XXO chain to solve it. Due to the simple form of the XXO basis, 
all further calculations can be done explicitly. The form factors of the model in 
the finite box are calculated in Section 2. The normalized mean values of the 
products of two canonical fields with respect to Bethe eigenstates are obtained in 
Section 3. In Section 4, the expressions for the two point temperature correlation 
functions in the thermodynamic limit are derived which is the main result of our 
paper. Important particular cases are discussed in Section 5. 
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1 The eigenstates of the Hamiltonian 



The model in the finite box of length L with the periodic boundary conditions 
is defined by the secondary quantized Hamiltonian 

L 

H = J dx{(d x ip + d x ilj) + c : (ip + *p) 2 : -h(^ + ^) + B(iP + a z ^)} , (1.1) 
o 

where c is a coupling constant, h is a chemical potential, B is a constant ex- 
ternal transverse field and : : means the normal ordering. The Pauli matrices 
are normalized as [a x ,cr y } = 2ia z . The one-dimensional fields ip a (x) and ipa( x ) 
(a = 1,2) are quantum operators in the Fock space with the canonical commu- 
tation relations for the Bose gas (e — +1 in the equations below) and with the 
canonical anticommutation relations for the Fermi gas (e = —1): 

^ a ( x )^p(y) ~ £ ^p(y)^c t ( x ) = S a pS(x - y), 

^(x)^(y)-e^(y)^+(x) = 0, ^ a (x)^(y) - e^(y)^ a (x) = 0. (1.2) 

In what follows we consider both cases simultaneously. It appears that the de- 
pendence of the correlation functions on the statistics is rather simple. 

The two-point temperature correlation functions are defined as the tempera- 
ture normalized mean values, 



Gl'\x,t;h,B) = — 



e-^+(x,t)^(0,0) 



Tr 



Tr [e- H / T ] 
a (x,f)^+(O,O) 



-HIT 



G^(x,t;h,B) = — t ";„.:. a ^ , (1.3) 



Tr [e~ H / T ] 
where 

tifrt) = e itH ^(x)e- itH , Mx,t) = e itH i; a (x)e- itH (1.4) 

and the traces are taken in the whole Fock space. Correlators G\ and G 2 are 
related by changing the sign of the magnetic field, 

G^\x,t;h,B) = G^\x,t;h,-B), (1.5) 

so that it is sufficient to calculate the correlators G^\x, t; h, B) only. Our aim is 
to calculate the correlators in the thermodynamic limit, i.e., at L — > oo keeping 
h and B fixed. To do this, we perform first the calculation for the finite box and 
then go to the limit. 

The basis in the Fock space of the model is constructed by acting with oper- 
ators i/j^i- 10 ) on t° t ne pseudovacuum |0) defined as 

ij a (x)\0) = 0, (0\ij+(x) = 0, (0|0) = 1. (1.6) 
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We say that a state belongs to the sector (JV, M) of the Fock space if it contains 
N — M particles of type 1 (a = 1) and M particles of type 2 (a = 2), so that 
the total number of particles is JV. The number of particles of each type being 
conserved separately, an eigenstate of the Hamiltonian can be obtained as a linear 
superposition of the basis states from the same sector. In the sector (JV, M) the 
eigenstates are enumerated by two sets, {k} = k±, . . . , k N and {A} = Ai, . . . , A^, 
of unequal (in each set separately) real numbers. Thus, the eigenstates in the 
sector (JV, M) can be written in the form 

L L 

\*N,M({kY, {A})) = J d Zl ...Jdz N Xn,m n (zu • • • , z N \{k}- {A}) 

o a\,...,a N =l,2 

x^i)---^ N M\0), (1.7) 

where the wave function xl*M tN ( z ii---i z N\{k}',{\}) is not equal to zero only 
if M elements are equal to 2 and N — M elements are equal to 1 in the set 
ai,..., a N , i.e., aj = N + M. 

The wave functions for the finite coupling constant c were obtained in ITD] (see 



also [jlql ). After imposing the periodic boundary conditions, the eigenf unctions 
of the Hamiltonian are given by the two-component (nested) Bethe Ansatz. The 
second component of the nested Ansatz (the "auxiliary lattice problem") is the 
Bethe Ansatz for the inhomogeneous XXX chain, which gives the higher weight 
eigenvectors. The complete basis is then obtained by acting with the Yangian 
generators onto the Bethe Ansatz vectors. 

The model at c = oo is defined by taking the limit c — > oo in the expressions 
for the wave functions. This can be done differently. One way to do this is 
to take the limit in the final expressions for the Bethe eigenf unctions, keeping 
the XXX chain eigenfunctions for auxiliary lattice problem. We use another 
way. It is to take the limit in the expressions for the wave functions before 
imposing the periodic boundary conditions, and to satisfy the periodic boundary 
conditions after that. On this way the auxiliary lattice problem is reduced to 
the eigenstate problem for the cyclic shift operator on the lattice which can be 
solved explicitly, e.g., using eigenfunctions of the XXO chain with the periodic 
boundary conditions. The advantage of this formulation is also the completeness 
of the XXO eigenfunctions in the isospin space. 

The Yang's wave functions for c = oo can be written in the form 

X a N \ M aN (z 1 ,...,z N \{k};{\}) 



1 

JV! 



E(-^) [P1 ^V PJV ({A}) e{z n <...< z PN ) 



det7v{e ifctl2b }, (1.8) 



where ^^^({A}) are components of a 2 iV -dimensional vector £at,m({A}). The 
sum in (|1.8|) is taken over all the permutations of X numbers, P : (1,2,..., X) — > 
(Pi, P2, . . . , P/v); [P] denotes the parity of the permutation. 
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The periodic boundary conditions for the wave function give the auxiliary 
lattice problem which is much simpler for c = oo than in the general case: 

QT^m) = (-e) N+1 e^ L C N 2 ; M aNai m)- (1.9) 

These equations should hold for all quasimomenta fci, . . . , fcjy. They can be re- 
garded as the eigenvalue problem for the cyclic shift operator Cn acting in a 
2 JV -dimensional vector space which can be identified with the space of states for 
the spin- 1 chain with N sites. 

The vectors £n,m, which have to be eigenvectors of cyclic shift operator Cn, 
can be chosen as the eigenvectors of the XXO Heisenberg spin-^ chain with the 
periodic boundary conditions. We write the expressions for £ N M in the form used 

m 



in 



N N 

Civ,M({A})= ■■■ 12 yjv,Af(ni,...,n A f|{A}) a^ ni) ■ ■ -% M )|ftiv), 

"1=1 «AI = 1 



<PN,M(ni,...,n M \{X}) = 



\{ sgn(n ; -n j ) 

l<j<l<N 



det M {e iXanb }, (1.10) 



where | f|\/y) = ®n=i\ In) ( an spins up) is the pseudovacuum for the chain 
with N sites, (ff-jy | ff-jy) = 1. The eigenvalues of Cn on the vectors (|1.10| ) are 
exp (i EtU A*) • 

The periodic boundary conditions for the function ifN,M together with equa- 
tion ( |1.9| ) result in the system of the Bethe equations, 



e ik a L = (_ e )iv-i e iEfci x \ a = 1, . . . , N, 



1.11) 



e i\ b N = b=l,...,M. 

The permitted values (k a )j G (— oo, +oo) of the quasimomenta k a (a = 1, . . . , N) 



are 

M 



2tt / A +e\ N - 1 \ 1 „ 



L V v 2 / 2 V L, 

0=1 

while the permitted values (Xb)i € (— 7T, 7r] of the quasimomenta A (6 = 1, . . . , M) 
are 

2n { N 1 + (-1)^-^ \ 

W^^f-y-^V — + / = 1 '---' iV - (i-i^) 

Each quasimomentum k a (|1.12|) is a function of the total quasimomentum Y,l=i ^6 
of the corresponding auxiliary lattice problem only. So one can decompose the 
quasimomentum k a into two parts 

1 M 

ka = ~ka+ T 12^b (1-14) 
^ 6=1 
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where k a does not depend on Aj, . . . , Am- We will often use this fact in what 
follows. 

The eigenvalues E^ t M({k}) of the Hamiltonian 



H \^N, M ({k}; {A})) = E NM ({k}) \y N M{k}; {A})), (1.15) 
are given as 

E NM ({k}) = £(^ 2 - ^ + B ) ~ 2MB - (1-16) 

3=1 

In the case considered (c = oo) they depend only on the total momentum of the 
auxiliary lattice problem which results in the additional degeneration of eigen- 
states in comparison with the case of finite c. In particular, the ground state at 
zero magnetic field B is degenerate. Indeed, the total momentum of the auxil- 
iary lattice problem for the ground state should be equal to zero at B = 0, and 
the requirement of the minimal energy fixes only the total number of particles 
N = Nx + N 2 , the numbers Ni = N — M and N 2 = M being otherwise arbitrary 
(of course, there exists also the additional degeneracy if even X and M are fixed). 
On the contrary, the ground state is not degenerate if B ^ 0. It is filled by the 
particles of the first kind {N\ = X, N 2 = 0) if B < and by the particles of the 
second kind (Xi = 0, N 2 = N) if B > 0. Thus, there is the phase transition at 
point B = 0. It should be mentioned that the quantization of the total momen- 
tum is the same for the XXX and XX0 chains so that the same results can be 
obtained using the XXX chain formulation of the auxiliary lattice problem. 

Consider now the properties of the eigenstates of the Hamiltonian. One can 
prove that two eigenstates belonging to different sectors as well as two eigenstates 
from the same sector but with different sets of quasimomenta are orthogonal, i.e., 

(VN'MWh W)\^N, M ({k}; {A})) = o (i.i7) 

if X' ^ X or M' ^ M or {k 1 } ^ {k} or {A'} ^ {A} (two sets are considered to 
be different if one of them cannot be obtained from the other by a permutation 
of the elements). The normalization of eigenstates is easily computed to be 

V&N,M{kl, • • • , k N ] Ai, . . . , \m)V& N,M\k\, . . . , k]y] Ai, . . . , Am)) 

= N M L N . (1.18) 
The eigenstates form the complete set, 

g£ ^ I^m(W;{A}))(^, m ({A:};{A})| _ ^ (Lig) 



N=0 M=0 fci<...<fcjv 
Xl<...<X M 



where the summations are performed over all different sets of solutions ki, . . . ,kjy 
( |TT2|) and Ai, . . . , A M (^JJ) of the Bethe equations ( |OT[ ). 
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Due to completeness of the eigenstates one can represent, e.g., the temperature 
correlation function G[~\x, t; h, B) in the form 

M-) (x t . hB) _ £e-^»(^(x,^(0,0)) jy , tf 

where (ipf (x, t)-0 1 (0, 0))n,m is the normalized mean value of the operator 
if>i(x, ^^(O, 0) with respect to the eigenstate \^N,M({k}; {A})) and the sum- 
mations in (|L2y) are performed as in (|1.19|) . The similar representation can be 
written for the correlator G± (x, t; h, B) as well. 



2 Form factors of local field operators. 

In this Section we calculate form factors of operators ipt(x, t) and i/;p(x,t), 
i.e., their matrix elements between two eigenstates of the Hamiltonian. Since the 
form factors of operators ipp(x,t) and ipp(x,t) are related by means of complex 
conjugation, it is sufficient to calculate the form factors of operators ipp(x,t). 
The nonvanishing form factors of if)p(x,t) are 

J%{x,t) = (^({g};{/i})|^(x,t)|^ +1 ,M({fc};{A})), (2.1) 
where the notation 

M=i M if/5 = 1 (22) 
M \ M - 1 if = 2 [Z - Z) 

is introduced. 

The form factor !F^ M {x,t) depends on the values of quasimomenta {k}, {A} 
and {q}, {fi} corresponding to the right and the left eigenstates in (j2.1|), respec- 
tively. Here and below the quasimomenta {k}, {A} correspond to the eigenstates 
in the sector (N + 1,M), and the quasimomenta {q}, {//} are prescribed to the 
eigenstates in the sector (JV, M). The Bethe equations for these momenta are 




e ik a L = U (- £ )N 


a = 


1, 


...,JV + 1 


e i\ b {N+l) _ (_l\M-l 


b = 


1- 


...,M 


e ig a L = V (_ £ }N-1 


0=1, 




,N 


e iV b N _ /-QAf-1 


6 = 1,. 




,M 



(2.3) 

where u and v are eigenvalues of the cyclic shift operators of the corresponding 
auxiliary problems: 

a; = exp^^A^, v = exp ( i J2 J • (2.4) 
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Using relation ( |2.5| ) and the eigenvalues ( |1.16| ) of the Hamiltonian, one ex- 
tracts the time dependence of the form factor (|2.1| ): 



!, N N+l s | 

^(X>? - E k " + hp) | ti&tW- ( 2 - 5 ) 

Here hi = h — B and h 2 = h + B are the chemical potentials of the particles of 
type 1 and of type 2, correspondingly. We denote 

T^ M {x)=T%{x,Q). (2.6) 

The distance dependence of the form factor is described by function ^^^(x). 
Using expressions (|1.7|) for eigenstates one represents it in terms of the wave 
functions as 

L L 

F ( N %(x) = (N + l)\ J d Zl . . . J dz N £ x a N 1 ^ aN (z 1 ,...,z N \{q};{ f i}) 





ai...a N /3 



oi,...,ojv 



x X a N 1 + 1 jn^---^N,x\{k};{X}). (2.7) 

Here and below the bar denotes the complex conjugation. Substituting now the 
explicit expressions for the wave functions into (|2.7|) and producing the necessary 
calculation one comes to the expression 



Fn?m( x ) = jfi J dz i ■ ■ ■ J dz N E \ e ( z Ri < ■■■< z Rn < x)F p 

RgSn \ 



(N) 



N-l 



+ E ( z Ri <---<z Rj <x< z Rj+1 <...< z Rn )(-s) n j F/3{j) 



+6(x <z Rl <...< z Rn )(-b) n F p (0)\ ]T (-l^e-^^-^N^) 

) Q&S N 

^2 (-i)\- p } e i{ - kp i Zl+ '" +kp N z N+ k p N+1 x ) ^ ^2.8) 



x 

P€S N+1 



where 

,.c*i. ..aj/3aj+i.. .a 

Ql,...,OjV 



ZTi f-\ — ^ ai ...Q N ( .a 1 ...a j f3a J+1 ...a N 



Further derivation requires more detailed explanation. 

First, note that £,n+i,m an d Cnm are eigenvectors of the cyclic shift operators 
Cn+i and Cn with eigenvalues ui and u, respectively (see ( |1.9|) and (|2.3|) , (|2.4j )), 
namely, 

j-a\a2-..a N+1 y.a2-..a N+1 ai ^aia2---ojv ..ta2-.-a.ivai 

S7V+1,M — ^ t,N+l,M ) s tv,M ~ ^N,M ' 

= 1, ^ = 1. (2.10) 
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Using relations ( |2.10| ) (N—j) times for both vectors, one can move the superscript 
(3 in ( p.9|) to the right, preserving the same order of indexes ax, . . . , of m 
and £,n+i,m with respect to which the summation in Q2.9Q is performed. This 
results in the expression 

Fp(j) = {G> V ) N -lF fi , Fp^Fp(N) (2.11) 



(one should take into account that u = u l ,u = v 1 ). Due to ( 2.11 ), one can 



move Fp in (|2.8| ) out of the braces leaving the factors {u)v) N ~i instead of Fp(j) 
in the first sum (over the permutations R G Sjy). 
Second, let us introduce the function 

p(z) = 9{z) - euv9{-z) (2.12) 

(which is an analogue of function sgn(z)). It is not difficult to see that the 
following relation holds 

{o{zri < ■ ■ ■ < z Rn < x) 

N-l 



+ °( z Ri < ■ ■ ■ < z Rj < x < z Rj+1 < ... < z Rn )(- 



eojv 



+9(x < z Rl < ...< z Rn )(-£uju) n ^ 



N 

= Hp(x- Zj ) (2.13) 
for zi, . . . ,zn,x G [0, L\ being all different. If two of z's do coincide, then the 



corresponding terms do not contribute into (|2.8 ) due to the fact that the last 
sum with respect to the permutations P G SV+i then vanishes. Hence the value 
p(0) is unessential. Due to ( p. 11] ) and (|2.13|) , the expression for M (x) can be 



rewritten in the form 

If h n 

Fn?m( x ) = -w F P J dz 1 ...J dz N H p(x - Zj ) 

n n 7=1 



X 



V" (_l^[PMQ] e i (kp 1 -QQ 1 )zi+---+i(kp N -qQ N )z N +ikp N+1 x ^ (2.14) 



Qes N 

P£S N+1 



Since the dependence on zi,...,z N in ( |2.14| ) is factorized, the integrals can be 
taken explicitly, 

( dzpix - z)e i{ - k - q)z = + eG)v)%, -. (2.15) 

{ (k-q) 
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The Bethe equations Q2.3| ) for k and q were used in derivation of ( |2.15| ). Due to 
this relation we have 



[pes N+1 K Pi 



4(l) = H)"(l :.^) A i^ V f-D r| — X-.-X 



Qi kp N - q N 



x 



expUx^k.-j^qM. (2.16) 



The expression ( [2.16|) is the desired representation for the form factor J r ^\ I (x). 
As usual, the dependence of the form factor on x is described by the "transla- 
tional" exponential. It is multiplied by the determinant of a matrix depending 
on the momenta {k} and {q} only (the sum in the first braces in ( |2.16j )). The 
factor Fp depends on the sets {/x} and {A} of the momenta of the auxiliary lattice 
problems only being a special "scalar product" of two Bethe vectors defined on 
the lattices with N and N + 1 sites: 



.,a N 



Now we are going to obtain the representation for Fp as a determinant of a 
matrix depending on {/z} and {A}. To do this it is suitable to use the notations 
for the vectors £ N m and £,n+i,m as the bra and ket vectors, (£jvmI an d |^jv+i,Af)j 
which are built upon the pseudovacua (f^Ar | := ®^ =1 (Tn | and | ft7v+i) := ®n=i\ In) 
of spin chains with N and N + 1 sites, respectively. Let us introduce the new 
vector 

{Zn,m\ =(£n,m\®(U+i\ (2.18) 

which is built upon the pseudovacuum (ft'jv+i | = ®^] 1 (Tn I- Then the scalar 
product ( J2.17|) can be rewritten as 

Fi = (6v,m(M)|W,m({A}), (2.19) 

F2 = [£,N,M-l({^}) °"(iV+l) %N+1,M 

({A})), (2.20) 

where the dependence on quasimomenta {/i} and {A} is written explicitly. Note 
that there are M quasimomenta in the set in ( |2.19| ) and M — l quasimomenta 
in the set {//} in ( |2.20| ); they satisfy different Bethe equations but nevertheless 



we denote these two sets by the same letter. At any particular case, it will be 
clear what the Bethe equations for the set {/i} are. 

As mentioned above, the factors F\ and F 2 given by ([2.19D and (|2.20 ) can be 



represented as determinants. This representation will be obtained now for F\. 
The derivation of the representation for F2 is quite similar. 
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Substituting the expressions for the Bethe vectors involved into ( |2.19| ) we get 



N N+l 



Fi = Y Y ^Ar,M(mi,...,mM|{A*})^Jv+i ) M(rii,---,nM|{A}) 

mi ,...,m,M =1 n\ ,...,nj^f=l 

x ( W \tf mi) • ■ • <r+ mM) (T ( - 1} ■ • ■ a ( - m) I W) • (2.21) 

Since there are no cr^v +1 ) here, all the terms containing &7 N+1 \ do not contribute 
in ( 2.21Q . Hence the summations with respect to n\ , . . . , um can be performed 
also from 1 to N. Thus (|2.21|) is rewritten as 

N 



F 1 = M\ Y ^N,M(n 1 ,...,n M \{fj})(pN+i,M('ri 1 ,...,n M \{X}) 



ni,...,n M =l 
1 N 

— JL ^_l'j[P] e -i[^P 1 ni+-+HP M n M ) 

ni ,...,n M =i PeS M 

x Y (-l) [Q ]e i(A<3ini+ '" +AQ M nM) , (2.22) 
which gives 

p x= (-1) [P] ( Y e l(Ap i^ tl)ni ] x ••• x [ Y e t{Xp M-^ n A . (2.23) 

PeS M \ni=l / \™ M =i / 

This is, by definition, the determinant of M x M-dimensional matrix. The de- 
terminant representation for F 2 can be derived in a similar way, only the terms 
with cr7 N+1 \ contributing. 

Let us now sum up the results of this Section. 

For the form factors J^ M {x^ we obtain the following representation 
■^LrOM) = (-*)"(! + eu>v) N det M B f3 det N+1 D 

fN+l N \ 

x exp I Y (~ itk l + ixk a ) - Y (~ it( l 2 a + ixQa) + ithp \ , (2.24) 

I a=l 0=1 J 

where hi = h — B, h% = h + B. The matrix elements of M x M-dimensional 
matrices B\^ are 

N 

OBiU = £e^-^>\ a,6=l,...,M, 

n=l 
V 

3 i(A„-/i t )ii 



(B 2 ) ai = Ee i(Vft) ", (S 2 ) 0i m = 1, 



n=l 



a = l,...,M, 6 = 1,...,M-1, (2.25) 
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and u, v are given by ( |2.4j ). The (N + 1) x (N + 1) -dimensional matrix D has 
matrix elements 



(D) ab = - , (D) aN+ i = 1, a = l,...,JV + l, 6 = 1,..., AT. (2.26) 

k a - q b 

It should be noted, that the determinant of the (N + 1) x (iV + 1) -dimensional 
matrix D can be represented in terms of the determinant of the NxN- dimensional 
matrix: 



detjv+i-D 



1 + — 

dz 



det7v(^i - zA 2 ) 



(2.27) 



2=0 

The N x iV-dimensional matrices A\ and A 2 have elements 



(A 1 ) ab = -^—, (A 2 ) ab = - , a,b = l,...,N. (2.28) 

k a - qb k N+1 - q b 

Since the matrix A 2 is of rank equal to one, detjv^i — zA 2 ) is a linear function 
of z. Thus, the r.h.s. of ( p.27| ) is equal to detA^Ai — A 2 ). The relation ( |2.27[ ) will 
be used in the next Section. 

Let us note that the permitted values of q a depend on the sum of the corre- 
sponding /i's. Hence the values of matrix elements of matrix D as well as the 
values of exponents in ( J2.24 ) are in fact different in cases (3 = 1 and (3 = 2 though 



being written down in the same way for the sake of simplicity. 

Let us consider the particular cases M = and M = N + 1. IfM = 
then J~^ M (x,t) vanishes due to its definition while J~^ M {x,t) coincides with the 
form factor of the one-component gas. In this case we put detAf=o-E>i = 1 (see 
the derivation of F\ above), v = u = 1, and the Bethe equations for {k} and 
{q} coincide with those in the one-component gas. Similarly, if M = N + 1 
then JN M (x,t) vanishes and ^F^ M (x,t) coincides with the form factor of the 
one-component gas (in this case Af, = 0, J2b=i /^-b = 0, hence v = uo = 1 and 

Bethe equations for {k} and {q} coincide with those in the one-component gas). 
The value of detM=N+iB 2 is a numerical factor related to the normalization of 
eigenstates (see ( |1.18|) ), being hence unessential. 



3 Normalized mean values of bilocal operators. 

In this section we derive the representations for the normalized mean values 

(^ + i,Af({fc};{A})|^(g,t)^(0,0)|^ +1 , M ({fc};{A})) 

(^ + i,M({fc};{A})|^ +liM ({fc};{A})) 1 • ' 
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and 



(^,Af(M;M)|^(x,t)^(Q,0)|^iv,M({g};{/x})) 

{*nA{&M)\*»A{&M)) 



(3.2) 



(P = 1, 2) with respect to the eigenstates of the Hamiltonian. These quantities 
enter the representations for temperature correlation functions. Below the deriva- 
tion of the representation in the case /3 = 1 is considered in detail, the derivation 
in the case (3 = 2 being quite similar. 

Consider the normalized mean value ( |3.1| ) for (3 = 1. Inserting the complete 
set of eigenstates (see Q1.19D ) between operators ipi(x,t) and ^(O, 0) and taking 
into account the normalization condition (fTTH) we get 



1 



(4>lfa ^^(O, 0))N+1,M - L 2N+l]\fM(]\f + l)M E Fn,m( X i 

' ll + ^l^ldetM.Bil^detyv+i.D) 2 



L 2N + 1 N M (N + 1) M 



fN+1 N \ 

x exp < ^2 (Hk* — ixk a ) — E(^& — ixqb) — ith\ \ . (3.3) 

I a=l 6=1 J 

Recall that u = exp iA and v = exp iQ where A and G are the total momenta of 
the corresponding auxiliary lattice problems: 

M M 

A = ^ A a , 9 = $> a . ( 3 - 4 ) 

a=l a=l 

The summation over {q}, {/j,} means the summation over all possible sets of quasi- 
momenta {q} and {fi}. Since each q a can be represented as q a = q a + Q/L (see 
( |1.14| ) and the comments there) where q a does not depend on /x's, one can perform 
the summation over {q} = q ly . . . ,q N independently of the summation over {//}. 
The expression under the sum in Q3.3| ) is symmetric under the permutations of 
g's (and g's as well) and /i's separately, being equal to zero whenever two g's or 
/i's coincide. Thus one can change the sum in (|3.3|) to the sum over all permitted 
values of each q a and each /x b , 

{q},{fi} 11<-«IN <?1 <?jV Ml MM 

M1<--'<MM 

where the sums over each individual q a (a = 1, . . . , N) and /i& (b = 1, ... , M) are 
independent: 

N 

E /(g.) = E /((ffa)i), E /to) = E /(«)■ (s.e) 
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The permitted values (q a )j and (nb)i are solutions of the Bethe equation 
Explicitly, 

t~\ 27r ( fl+e\N-l A . „ 

(fc)i = T (-(— ) — + J J' JGZ ' 
2n f N 1 + f-D^-A^ \ 

W = f (-y-^V — + l = l --- N - ( 3 - 7 ) 

Taking into account ( |3.5| ), let us perform first the summations over q\, . . . , 
and then over jjli, . . . , [Im hi ( |3.3| ). Note that the dependence on g's in ( |3.3| ) enters 
only the elements of matrix D and the exponential factor. The determinant of 
the matrix D can be represented as the sum over permutations, 

det N+1 D= £ (-1)11 J] (3.8) 

Pg5jv=i a=l K Pa Qa 



so that the dependence on each q a in the r.h.s. of ( j3.3|) is factorized. Therefore we 
can use the technique similar to that applied in the cases of the one-component 
impenetrable Bose-gas (16], [I7j and of XXO model [ 18), . This procedure (which 
can be called "inserting the summation into determinant" ) results in our case in 
the following indentity: 

1 1 + evuj 



'U}\ 2N ( N+1 N 1 

qi (jjv v a= l 6=1 ' 

(3.9) 



L 2N+ 

qi <?jv 



2=0 

The matrix elements of the (N + 1) x (JV + 1) matrices and -R^ - -* are 

M _|_ £j/ q|2 e -itq 2 +ixq ^ Q 

(it^U = e - (fc ^- (fcfe) , (3.10) 
the function e-(k a ) is defined as 

e_(fc a )=expf^^ -J (3.11) 

and 2; is a parameter. Since the rank of matrix R(~) is equal to one, det tv+i (S^ + 
zR(~}) is a linear function of z. 

To go to the thermodynamic limit (which is done in the next Section) it is 
necessary to rewrite the expressions for the matrix elements of the matrix 
Note that matrix elements given by ( |3.10|) for finite L are well-defined functions, 
since all possible zeros of the denominator of the expression under the sum over q 
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are canceled by zeros of the numerator coming from the factor |l+ez/u)| 2 . It is due 
to the fact that the condition k — q can be satisfied only if A — 6 = 7rii £ (mod27r), 
which follows from the Bethe equations. In the limit L — > oo, however, the values 
of k and q become arbitrary, the summation over q is changed for the integration 
j J2g — » 2^ I^°oo dq and poles on the integration contour appear which should be 
taken into account before going to the limit L — > oo. Therefore one has to rewrite 
the matrix elements of in terms of functions which are well defined in the 
thermodynamic limit. 

The necessary calculations are given in Appendix. As the result, the matrix 
S (-) can be represented in the form 

5 H = J + l + ^os(A-0) Vt ] + " MA 2 ~ 9) Vt\ (3.12) 
where I is the unit (N + 1) x (JV + 1) matrix, (J) afe = 5 a6 , and matrices V} 2 are 
, T/( -K 2 ef } (A; a )e_(A;&) - e_(A; a )ei" ) (A; 6 ) 

U "Z • (3 ' 13) 

The functions involved are defined as 

e H(jfe) = eH(jfe) e _(jfe) > e _(jfe) = e ^%^ (3 14 ) 

o —itq 2 +ixq —itk 2 +ixk r\ 

e { -\k) = jY.- 1 . <7 = <7 + y- (3-15) 

These functions are well defined in the thermodynamic limit. 

The summation over g's being fulfilled ( |3.9|) , one remains with the sum over 
the quasimomenta /x's in the expression (|3.3|) for the normalized mean value: 



(1>t{x, t)Vi(0, 0)) N+ 1,M = ArM(Ar+1)MM , E • • • E l det A/ 5 ll 



(3.16) 



2=0 



Let us perform the summation over . . . ,^m- Note that each element of 
matrix S^' depends on the sum 6 = (ii + • • • + (jlm only. Moreover, since 
the summation over q in ( 3.1U| ) is governed by the rules ( |3.6| ), ( p.7| ) and since 
q = q + Q/L, matrix elements (S) a b are periodic functions of 6 with the period 
2n. Matrix does not depend on fii,...,fi M at all, hence detjv+i^ - ^ + 
zR(~') is also a periodic function of 0. Further, expz© = v (see ( |2.4j )) is an 
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eigenvalue of cyclic shift operator Cn, hence G takes the values 27m /N where 
n = 0,1, . . . , N — 1 (mod N) and we can write 



N-X 

1 = E 6 (N) 

n=0 



N 



fix + • • • + Hm 
2^ 



n 



where 8(n) is Kronecker symbol on Zn, defined as 
m E Z, 6(N)(m) = 



1 if to = (mod N) 
otherwise 



and can be represented as Fourier sum 



-i N-X 

1 x — T 2ttz 



'(TV) 



TO 



iV 



E 

p=0 



e w 



(3.17) 



(3.18) 



(3.19) 



Inserting ( |3.17| ) under the sum over fix, . . . ,\im in ( |3.16|) we can change — > 
27m/ TV everywhere in the expressions for elements of the matrix S^~\ Using the 
representation ( |3.19| ) we can rewrite ( |3.16| ) as 



(V> 1 + (M)V'i(o,o)W,A f 
l 



jV-1 



N M + 1 (N +l) M M\ pt 
d 



E-EE 



Jp(fll-\ \-HM)-^rPn 



|det M -Bi 



xe 



-ith-i 



dz 



Mm n,p=0 

det N+1 {S^ + zRV) 



(3.20) 



2=0 



where matrices ) now do not depend on fix, ■ ■ ■ , fJ>M at ah and their elements 
are (S}~') a b = (S^) a b\e=2nn/N- Matrices have the structure corresponding 

to CT) , 



I + £ cos(A - ^ 



£ sin(A - 2h> 



(-) 



2 5 



(3.21) 



2 i,n 2 

|e=27rn/Ar, i-e., one should put q = q + jtt in the function 



where V$ = V± 

(|3.15[) determining the elements of the matrix V - / . 

Considering the determinant of matrix B\ as the sum over permutations, 



M 

det^ = x: (-i) [p] n 

PeS M a=l 



N 

E • 

n a =l 



(3.22) 



we see that the dependence on fix, ■ ■ ■ , fiM in ( p.20| ) can be factorized. Thus 
we can perform the procedure of "inserting summation into determinant" with 
respect to the summations over fix, ■ ■ ■ , fiM- Picking up the relevant terms from 
( |3.20| ) one comes to the indentity 

1 



N M (N + 1) M M\ ^ 



. . . e «p(MiH — i-mat) 



det a/-Bi I 



detMf/i 1, ] : 



(3.23) 
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where M x M matrices U^ 1 ' ' have elements 

i N N 

(f/ (i,-) U = E E E e^ + ™-^ +inX °- imX K (3.24) 

iV -+- L) n=1 m=1 ^ 

Finally, the normalized mean value of operator ipf (x, t)?/>i(0, 0) takes the fol- 
lowing form 

i N-l 

(^ 1 + (a;,^ 1 (0,0)) J v +1 ,M = e-^- E ^ pn 

iV n,p=0 

xdetuU^ [det N+1 (S^ + _ det^^H] , (3.25) 

where the property Jj/(z)| z= o = /(l) — /(0) valid for a linear function /(z) is 
used. Expression ( |3.25| ) is the final expression for the normalized mean value 
considered. 

For the normalized mean value (^(x, t)ip 2 (0, 0))n+i,m au the calculations 
presented above can be done similarly. The main difference is in the summation 
over the values of quasimomenta of the auxiliary lattice problem, i.e., over the 
set {/i}, which in the latter case is {//} = //i, . . . , \im-\- The representation 
analogous to equation (|3.25|) for this normalized mean value is 

i N-l 

(^( x ,t)^ 2 (0,0)) N+1M = e-^- E e"^" 

^ V n,p=0 

x 



det M (U^+P^)-det M U^\ 

det N+1 (S^ + flH) - detjv+^H] . (3.26) 
Here the M x M matrix pv~) is of rank one 

(**-')- = (3 - 27) 

and matrices Ujf'^ have the elements 

i at v 

(°? , " ) )«* = W^TTT E EE el(p+m_n (3.28) 
iv {IV -+- i) n=1 m=1 p 

It should be noted that though the matrices U^ 2 ~^ and U^~' look formally the 
same, their matrix elements are different (except, e.g., the case p = 0, see below) 
since quasimomenta {//} in these two cases satisfy different Bethe equations. 

Consider the particular = 0, t = 0. Then S n = I and only the terms 

with p = contribute into the normalized mean values. It is easy to see that 

(U£'- ] U = (U^U = 6 ab - (3.29) 
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and it follows from (Q), (Q) that 

(N + l) -M 



(iP+(0, 0)^(0, 0)) n+1 ,m 
(^+(0,0)^(0,0)^+1^ 



.1/ 



(3.30) 



which is quite obvious result since in this case the normalized mean values are 
just the expectation values of density operators for particles of type 1 and type 2, 
respectively. 

Let us consider now another normalized mean value, given by Q3.2|) . Inserting 
the complete set of eigenstates we get 



L 2N+i N M^ N + iyi ^E^ F$ m (x, t) F$ M (0, 0) 



L 2N + l N M (N + 1) M 



\l + evQ\ 2N \&A M B l 

{fe},{A} 



xdetAr+i-D 



1 + 



d_ 

dz 



detAr(v4i - zA 2 ] 



2=0 



X 



( N N+l \ 

exp <^ Y( it( lb - ix 1b) - E ( itk l ~ ixk a) + ithl \ ' 
I 6=1 a=l J 



(3.31) 



where the relation (|2.27|) is used for one of the form factors. Again, as above, we 
can perform the summation as follows 



E 

{k},{\} 



E 



fcl<--<fcjV+i ("^ ' ;. 

X-L<-<X M 



1 1 



(3.32) 



The procedure of "inserting the summation under determinant" with respect to 
the summation over ki, . . . , kw+i gives the following result: 



\l + evuj\ 2N 



#w+i(jV + l)! 4 



E " " ' E detjv+i-D 



ki 



1 + 



d_ 

dz 



detjv(y4i — zA 2 



2=0 



X 



( N N+l 

exp <^ E(^?6 - - E (^ fc a _ ^« 

I 6=1 a=l 



g(x,t) + 



dz 



det N (S {+) - zR (+) ) 



2=0 



Here we introduce the function 

k 



k = k+y. 



(3.33) 



(3.34) 
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The N x N matrices S^ + ' and are 

|1 +eisoj\ 2 



(S«U = e_(g a ) e_(g 6 ) 
|1 + evoo\ 2 



L 2 



E 



-itk 2 +ixk 



(lb 



L 3 



-(g.)E 



—itk 2 +ixk 

k-q a 



(k - q a )(k - q b ) 



k = k + y, 



-itk' 2 +ixk' 



. (3.35) 



The matrix is of rank equal to one, and detAr(5' < ' + ' ) — zR^) is a linear 
function of the parameter z. 

These matrices can be put into the form (see Appendix) 



S« = 1 + 



1 + e cos(A - 0) 



(+) 



e sin(A - 0) 



V, 



(+) 



2 5 



R (+) = 1 + e cos(A - 0) 



-s sin(A - 0) 4 +) + 1 = £ C °: (A = 9) R {+) 



3 5 



where matrices involved are 

/ T/ (+h 2 ei +) (g a )e_(g 6 ) -e„(g a )ei +) (g 6 ) 

1^1 Jab 



(W +) )a6 



(4 +) )ab 



(-R3 +) )a6 



£ <?a - g& 

2 [e_(g a )]- 1 e_(g h )-e_(g a )[e_(g & )]- 1 

L q a -Qb 

e ( :\qa)e { :\q b ) 



1 

2L 



e { :\q a ) | ej +) (g fe ) 
e-(gt) e_(g a ) 
1 



Le_(g a )e_(g b )" 



Here 



and 



e + +) (<?) = e-(g)e + (g), e_(g) = e*" 2 



■itq —ixq 



,(+) 



-itk 2 +ixk 



„—itq z +ixq _ A 

■ *=* + r 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



Due to ( |3.33| ), the normalized mean value (|3.31|) can be rewritten as 

1 _ _ . _ l2 



(0,0))^ 



iv^(iv + i)^5"'g |detA/jBl1 



xe 



det N {S (+) - 2i? (+) ) 



(3.40) 



z=0 
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The procedure of "inserting the summation into determinant" with respect to A's 
is performed in a similar way as above. Using the fact that detjv+iOS' — ^-R ) is 
a In periodic function of A and that A = |^ where m — 0, 1, . . . , N (mod N+ 1) 
we insert 

N 



1=-£5 {N+I) (m-(N + 1) 



m=0 



Ai + • — h Am 
2?r 



2nm 
N+l 



(3.41) 



in the elements of matrices 



under the sums over A's in ( |3.40| ). Changing A 
S( + \ R( + ^ and in function g(x,t) (|3.34j ) and using the representation ( |3.19| ) for 
the Kronecker symbol, we rewrite (|3.40|) as 

(^(x, 0^(0,0))^ 
1 



N 



N M (N + 1) M+1 M! ~^ f^r^o l ae ^i| 



xe 



ith\ 



g m {x,t) + 



d_ 

dz 



(3.42) 



2=0 



where subscript m means that everywhere in the functions determining the el- 
ements of the matrices S^ + \ and in the function g(x,t) one should put 
2Tcm/(N + 1) instead of A. 

Similarly to fl3.23j ), we have the relation 

iV M ( ^| 1) M M! E---Ee- fr(Al+ "- +A " ) Idet^l 2 = det M ^ +) , (3.43) 



where the M x M matrices U^ 1 '^ have elements 

i N N 



W 1 



+)' 



ab 



T) EEE 



-i(r+m— n) A+in/i a —im^ 



(3.44) 



Finally, for the normalized mean value (ipi(x, t)ipf (0, 0))jv,A.r we obtain the 
representation 

1 AT _ . 



(^ 1 (*,t)V'i + (0,0)> WlW = e^— 7 E ^ 



r,m=0 



X 



det M t/ r (1 ' +) fdet^W - + (g m (x,t) - 1) det iV ^ +) 



(3.45) 



For the sake of completeness, let us write down the representation for nor- 
malized mean value of operator ^ 2 (x, t)^(0, 0): 

i AT _ . 



(ifi 2 (x,t)ifi£(0,0)) NtM -i = e 



ith\ 



N + 



T E 



eiv+i' 



r,m=0 



X 

X 



det M _ 1 (C/ r ( 2 '+) - P r W) + (5 r , - 1) detM.!^ 2 '^ 
det JV (5'i + ) - i?L +) ) + - 1) detjv^W 



(3.46) 
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Here the (M — 1) x (M — 1) matrices U^ 2 ' + ' have elements (recall that in this 
case there are M — 1 /i's in the set {/i}) 



i N N 

( U ?' +) )<* = mN , u E E E e -« r+m - n » +im ^- m ^ (3.47) 
and the matrix Pj + - ) is of rank equal to one 

i N N 

{P ( + ) U = Yl E e-**-** x -rt £ E e-^'+^'-^. (3.48) 

iV {i\ tij A m=1 A , n=1 

The results obtained in this Section allow us to calculate the correlation func- 
tions in the thermodynamic limit. 



4 Temperature correlators in the thermodynamic limit 

Here we present the derivation of the temperature correlation functions 
( |1.3|) in the thermodynamic limit L —>■ oo. The correlation functions and 
(j?2 are related by inversing the sign of the external field, 

G^\x,t;h,B) = G [ ^\x,t;h,-B), (4.1) 

so that it is sufficient to calculate only one of them (an independent calculation 
of both correlators confirm, of course, the relation ( j4.1| )). Our result is the deter- 
minant representation of the correlators which is the generalization of the results 



obtained earlier for one-component models [16, 17, 18, 19 



Consider first the partition function of the gas which is the denominator in 
11.31). It can be written in the form 



oo N 

Tr' 



e -H/T 



1 + E E E E e-^-« 9 })/T (4 2) 

JV=1 M=0 qi<...<q N W<-<MM 



where the energy in the sector (N, M) is given by 



N 

EnMU}) =Y,(<l 2 a- h + B )- 2MB - ( 4 -3) 

a=l 

As usual, 

q a = q a + j^, @ = Mi H h Hm- (4.4) 

For finite L, the partition function can be presented as an explicit linear combi- 
nation of several determinants of "discrete" infinite dimensional matrices. This 
representation allows us to obtain strictly the thermodynamic limit [L — > oo) of 
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the partition function and also calculate the 1/L corrections to the free energy. 
The corresponding results will be given in a separate publication. Below we will 
consider the thermodynamic limit only. 

Note that due to the summation over all possible values of q\ < . . . < in 
( [4.2| ) it is possible, by making the same shift of all q a , q a — > q' a , to present each 
q a as 



Qa = q a + J, 



lei <_ I 



(4.5) 



So in the thermodynamic limit one can neglect in (|4.4j ). For the same reason 
one can neglect in the limit also the difference between the permitted values of 
q a in the sectors with odd or even number of particles N (see ( [1.12j )). Under this 
conditions the summation over fix < . . . < fiM and over M can be done with the 
result 



TV 



E E 

M=0 fll<--<HM 



;i + e-) 



The permitted values of quasimomenta q a (a = 1, . . . , N), 

T' 



(q a )j+x - {q a ) 



(4.6) 



(4.7) 



fill the interval (— oo, oo) densely in the thermodynamic limit and the sum with 
respect to q a should be changed for the integral 



7E/(9) = rE/((9)i) 



1 

2^ 



dqf{q). 



Thus, we have for the partition function 



Tr [e- H ' T ] = 1 + E E (1 + e^) N e-^^- h+B ^ 

N=Xq 1 <...<q N 



n 



1 + 2 cosh f e' 3 ^ 



(4.8) 



exp 



L 
2^ 



dq In I 1 + 2 cosh ^ e t 



L 



oo. 



(4.9) 



Our purpose in what follows is to present the correlators as Fredholm determi- 
nants of linear integral operators. The typical expression for both the numerator 
and the denominator (the partition function) in ( |1.3|) is of the form 



F=E E det N A=J2 £ 

N=0 <?i<.. .<<7iv N=0 qi<...<QN 



A{q x ,qx) 
A(lN,qx) 



A(q N ,q N ) 



(4.10) 
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(by definition, deto^4 = 1 for any matrix A). The sum in ( [4.10| ) is taken over 



all permitted values of ?i,...,9jv- m the thermodynamic limit, as explained 
above, one can neglect the difference between the permitted values of g's in the 
sector (N,M) and put, e.g., (q)j = -j^-j (in other words, it is essential only that 
(q)j+i — (q)j = x)- Then, due to (|4.8|), in the limit L — > oo the expansion ( |4.1U| ) 
is just the expansion for the Fredholm determinant, 

oo - 

F = H / dq 1 - ■ -dq N det N A = det{I + A). (4.11) 

N=0 q 1 <...<q N 

Here, by definition, we denote det ^A the determinants of the N x iV-dimensional 
matrices with matrix elements A(q a , qb)'- 



det nA : 



A(qi,qi) ■•• A(q u q N ) 



(4.12) 



A{q N , qi) ■■■ A{q N , q N ) 
Functions A(q a , qb) are obtained from the elements A(q a ,q b ) as 

A{q a ,q b ) = Km -^A(q a ,q b ). (4.13) 

The linear integral operator A acts on functions f(q) as 

oo 

(Af)(q)= J dq'A(q,q')f(q') (4.14) 

— oo 

and the kernel A(q, q') is given just by the equation ( |4.13|) . 

It should be kept in mind that in our notations A(q a ,q b ) denote matrix ele- 
ments of matrices A entering the expansions of the kind ( [4.1 0| ) which contain the 
sums. Quantities A(q a , qb) denote matrix elements of matrices entering the ex- 
pressions similar to ( |4.11| ) containing the integrals. Finally, A denote the integral 
operator with the kernel A(q, q'). 

The Fredholm determinant is well defined if the trace of operator A is finite 
/ dqA(q, q) < oo. It is the case for all the operators considered below with the 
exception of the integral operator corresponding to the partition function which 
is divergent in the limit (see ( |4.9| )). The partition function is the denominator 
in the expressions (|1.3| ) for the correlation functions. The mean values in the 
numerators are also divergent. It will be seen that this divergency is described 
by the same operator as in the partition function, so that the final answer for the 
correlator is finite. To make sense to the intermediate formulae, one should keep 
this divergent operator regularized. We will make this regularization writing for 
the partition function at L — > oo 



Tr 



e -H/T 



det(J + Z), Z{q,q') = 2 cosh f e-^5 L (q-q'), (4.15) 
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Here S^q — q') is a regularization of the Dirac delta-function; e.g., one can put 
sin L(q-q') 

Of course, this regularization reproduces only the leading term of the partition 
function; the 1/L corrections to the free energy should be calculated from the 
exact expression (|4.2|) . 

Turn now to the correlator G[ ^ (|1.3| ). Let us calculate the numerator of 
G{~ ] . There are three main steps in the calculation: i) take the thermodynamic 
limit; ii) sum over A's (or /i's) and M; iii) sum over N by means of the Fredholm 
determinant formula ( |4.11|) , extracting explicitly the determinant det(J + Z) 
(which cancel exactly the denominator). The result after taking L — > oo is the 
well defined Fredholm determinant representation for the correlation function. 

Consider the numerator of G[ . Using the representation ( 3.25 ) for the nor- 
malized mean value involved, we get 

Tr [e~ H/T ^f (a;,t)Vi(0,0) 

°° N+1 E N+1M ({k}) 

= EE E E c ' T — (Vi + (^,*)V'i(o,o)) J v + i ) m 

N=0 M=0 k!<...<k N+1 Ai<...<A M 

00 1 N-l N+1 

= e- U{h ~ B) E jj E ^ m E ^ M E det M U^ 

N=0 iV p,n=0 M=0 \ 1 <...<X M 

x [det N+1 (S^+R^)-det N+1 S^], (4.17) 

fel<...<fejv+l 

where tildes over matrices mean that the temperature factors exp{— Y,a=i(^a ~ 
h + B)/T} are included into the determinants, i.e., the matrix elements ( ) a b of 
the corresponding matrix is multiplied by the factor exp{ — ka ~!^ +B }. It means 
that the matrix §{~^ has the structure 



11 



kJ-h+B 



(SPU = e-^^5 ab + ■■■, a,b=l,...,N + l. (4.18) 

Since in ( [4.17|) the summation is over all values k\ < . . . < k^+i, there is the 
periodicity in A (recall that k a = k a + A/L) and one can insert the Kronecker 
symbol on Zn+i ( |3.41|) . 

Then in the thermodynamic limit, L — > 00, the numerator of G± acquire the 
form 

e-^ T ^ + (x,t)^(0,0)^ 

00 1 N N-l 

= e -u(h-B) \p 1 ST e^Ti rm -^r pn 
n=o N(N + 1) r m=0p n=0 

xX N>p>r J dk x - ■ -dkN+iEfcl^rfan), (4.19) 

fci<...<fejv+i 
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where we denote 

N+l 

X N , P ,r ■= E e¥M E det M f^;-\ 

M=0 Ai<...<A M 

2fc"+i(^n, m ) := detjv+i^^) - det^S^^). (4.20) 

The matrices (U,^.~^) have matrix elements 

= e^M^U (4.21) 

and elements of matrices (t^ 1 ' - -*) are given by ( p. 24 ). The matrix S^~\rj nim ) has 
the structure 

o { '{r]n,m) = N H Vi H V a> (4.22) 

where 

27im 2nn 

Vn,m := ]7TI " -j^- (4-23) 
The matrix elements of the matrix M are 

AA(fc a , fc 6 ) = e — 5 L (A; a - fc 6 ) (4.24) 



and the elements of the matrices Vip, > are obtained in accordance with the 
rule (|4T4]) : 

« . v E+(ka)E4k b ) - E.(k a )E + (k b ) 

Vl{Ka, Kb) — C T 



V 2 (k a ,k b 



n(k a - k b ) 

n_±B [E„{k a )]- l E„(k h ) - E^iE-ih)}- 1 



TV 



Ki-\KM)=e-^ E - {k i E - {h) . (4.25) 

(Since the elements of the matrices v[J and V12 , which are obtained from 
(Vx 2 )ab an d (^12 )o6 by means of (|4.14|) , are described by the same functions, 
there is no need to distinguish between Vj-^ and V^, hence we write simply 
Vi,2-) The functions E± are 

E + (k) = E(k) E_(k), E_(k) = e^, (4.26) 

where 

°9 p —itq 2 +ixq 
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These functions are the thermodynamic limits of the corresponding functions 
e+ , e_, which is explained in detail in Appendix. 



Consider the quantity X NtP>r . Note that, by definition, det^C/^ ) 



and it is easy to see that det ) 



M=0 



M=N+l 



0. Thus 



N 



X 



N,p,r 



1+E^ M E det M ^;-). 

M=l Ai<...<A M 



(4.28) 



To obtain a useful expression for Xn, p ,t let us note that the following relation 
is valid 



J2 det M U^ 

Xi<...<X M 



E 

l<mi<...<m M <N 



p,r 5 



(4.29) 



where 



detMM/W = E ("1) [P1 
Pes M 



M 
a=l 



N 



n„=l 



1 



? n VJV+1 ? 



7 -i\(r+m Pa -n a ) 



1 

AT 



e i^(p+m a -n a ) 



(4.30) 



To see that the relation ([4.29 ) holds it is sufficient to note that: i) det mU^~' 
is symmetric with respect to permutations of A's being equal to zero whenever 
two of them coincide; ii) detj^W^ 1 ) is symmetric with respect to permutations of 
m's (being zero whenever two of them coincide). After changing J2\ 1 < <a to 
1st ' ' ' i n the l.h.s. of ( f4.29| ) (and similarly for the summation over m's 
in the r.h.s.) the relation ( 4.29|) becomes obvious. 

Let us find the explicit expression for the elements of the matrix W^} . Using 
Bethe equations for A's and /x's we have the following summation formulae (here 
m := mi, . . . , m M ; n := m, . . . , n M ): 



N + 



3 — iA(r+m— n) 



r+m— 



(r + 7TI < N + 1) + 



\M+1 



r+m—n—N- 



_i#(r + m > iV + 1), 



_ \ " p in(p+m-n) 



5 p+m - n 6(p + m<N) + (-l) M+1 5 p+m ^ N 6(p + m > N) 



(4.31) 



These expressions are valid only for the values of the integers involved: p = 
0, . . . , N — 1; r — 0, . . . , N; mi, ... , mjw, ni, . . . , % = 1, . . . , N. Here and below 
we use the short notation for the usual Kronecker symbol, S a = S a> o. Making use 
of ( |4.31|) and summing over m, . . . , un in ( |4.30| ) we get 



J p—r+m a -m b 



1 



,M+1; 



Jp-r+ma— m b +N+l 



d(p + m a <N) 



4" |^<5jd— r+m. a - 



m b +l t \ i-J 0p- r + ma — mb -]\[ 
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9(p + m a > N). 



(4.32) 



Using the indentity for the determinant of an N x N matrix, 

H + det N A, 



N N 

det JV (/ + A) = l + J2 A Jd+ E 



3=1 



ii,i2=l 

J1<J2 



4 • ■ A ■ ■ 

^-32,31 ^J2J2 



(4.33) 



we obtain, due to the relation ( |4.29| ), that 

1 
2 



X 



N,p,r 



det N (I + e~W_) + det N (I - e~W. 
det N (I + e^W + ) - det N (I - e^~W + ) 
where N x N matrices W± are 



1 

+ 2 



(4.34) 



(W±) ab = I5p- r+a - b ± 6 p - r+a - b+N+1 ] 9(p + a< N) 

+ [5 p _ r+a _ 6+ i ± <5p- r+a -6-7v] <?(P + a > N). 



(4.35) 



To proceed further, let us look at the quantity z-N+i(i]n,m) as a function of 
r] ntm . Obviously, 

N 

Z { N Uvn, m ) = E e^-$H, (4.36) 

s=-iV 

where are some combinations of the elements of the matrices J\f, Vi,2,7^~' ) 
which do not depend on n, m. After summation over n and m in ( [4.19| ) the 
terms with s = 0, . . . , N in ( f4.36|) will produce <5 p _ r while the terms with s = 
—N, . . . , — 1 will produce <5 p _ r+ i. Hence we should take into account only the 
terms in ( f4.34|) for which either the condition r = p or the condition r = p + 1 is 
satisfied. It follows from ( f4.34j ) and ( f4.35p that 



X 



N,p,r 



r=p 



N,p,r\ r=p+ i 



(4.37) 



(4.38) 



1 + e 2B / T ) N -p y X~ I - fu- *2B/2V 
Let us introduce the notation 

7 = 1 + e 2B ' T . 
Then we can write the answer 

Xn, p ,t = 1 + (~f N ~ p — l)5p_ r + (7 P — l)5 P -r+i + unessential terms (4.39) 

where "unessential terms" are terms constructed from the Kronecker symbols 
containing in their arguments p, r and which are different from those written in 
( |4.39|) explicitly. These "unessential terms" will not contribute after the summa- 
tion over n, m,p, r in (|4.19|) . 
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Substituting ( |4.39| ) into Q4.19| ) we obtain the following expression under the 
sum over N and (N + l)-fold integration over fc's (except the overall factor 
exp{-it(h - B)}) 



1 



N N-l 

(^+i)^J?o£i 7 

N-l 

+ E(7 7V - p -i; 
P =i 



N 



^■Kir\ n ,raV _|_ g— 27TIJJ 



>s! v j 1 (r /n , m )+S« 1 (0). (4.40) 



Due to the explicit dependence on N of the elements of the matrix under the 
determinant in ( |4.19| ), some efforts are needed in order to apply the Fredholm 
determinant formula ( [4.1 1| ) to the expression Q4.19Q as a whole. The first trick is 
to change the double sum (over n, m) in the first term in (|4.40|) to the ordinary 
integral, i.e., the following indentity holds 



N N-l 



N-l 

+ E(7 iV - p -i; 
P =i 



+ e 



-2iririn, m p 



1=1 N +l{Vn,m) 



N 



N-l 

i + j2 (i N ~ p 
P =i 



e i«P + e -im> ( 4 - 41 ) 



It is sufficient to sum over n,m in the l.h.s. of Q4.41 ) and to integrate over r\ in 
the r.h.s. in order to see that the results coincide. 

The second trick is to expand the summation with respect to p in the r.h.s. 
of ( |4.41| ) to the infinity. Adding the term at p = N obviously do not produce any 
additional terms to the r.h.s. of (|4.41|) because ( r y N ~ p — 1)\ p =n = 0. The terms 
with p > N will die after taking the integration with respect to 77 since then the 
quantity e^Hjv^) (or e~ ir,p, E]sf(r])) contains e lv only in integer positive (or only 
in integer negative) powers, due to ( |4.3(j| ). 

Therefore, the expression (|4.40|) can be rewritten as 



1 



" . ! + J2(j N ~ p - 1) 
P =i 



:,"IP 



+ e 



-11)P 



=■(-) 



+ e 



-ir)p 



:,i'IP 



+ e 



-irjp 



^nIi(t,v)- 



(4.42) 



P =i 



Here we denote by H^ 1 (7; rj) the expression for E^l^r]) given by ( |4.20| ) where 
the replacement S^~\rj) — > jS^~\rj) has been made. The last expression in 
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(|4.42|) is what we need in order to apply the Fredholm formula ( |4.11| ). Taking 
into account that jj\f(k a , k b ) = Z(k a , kb), we get 



Tr 



-H/T. 



^ + (x,t)^(0,0)l =e^ h - B ^— / drjF^rj) 



2tt 



x 



det [I + Z + jV{r)) ) - det ( I + Z + -yV{r]) 



(4.43) 



Here 



~, . 1 + e cosr] ~ sinr/ ~ 



2 2 

and Z, Vi,2,7£ < ~' ) are integral operators with the kernels given by (|4.15|) , (K 
We have introduced the function 



(4.44) 



F( 7 , 77) = 1 + l~ p (e iVP + e~ lvp ) 
P =i 



(4.45) 



It is worth mentioning that 1 < 7 < 00 for any real external field B. At the 
point 7 = 1 (B = —00) the function F( 7 , 77) is equal (up to the factor 2n) to the 
27r-per iodic delta- function of the variable 77. 

Finally, taking into account the representation ( 4.15 ) for Tr[e~^/ T ] we obtain 
the following representation for the temperature correlation function considered 



det (/ + 7 V{rj) + H^) - det (/ + 7 V( V ) 



where integral operator V(r/) is of the form 
v 1 + e cosri - sinr/ ~ 

The integral operators Vi j2 , do not depend on 77 and possess kernels 
El(k)E^(k') - El{k')El{k) 



(4.46) 



(4.47) 



Vi(fc, fc') 
V 2 (A;, 



7r(fc - jfc') 



7r(fc - k') 

^. E T {k , } _ E T (k) m 



El(k) 



E T _{k' 



27T 

The functions E+(k), E^(k) are 



(4.48) 



(4.49) 
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and the function E(k) is given by ( |4.27| ). The Fermi weight $(k) is 



•&(k) 



-B/T 



2 cosh j; + e t 



(4.50) 



Let us now consider another temperature correlation function Gi ( |1.3|) . Us- 
ing the representation ( |3.45| ) for the normalized mean value, we represent the 
numerator of the temperature correlation function as follows 

00 N % M (W) 

= EE E E e — ^(^(x,t)^ + (o,o))^M 

N=0 M=0 q-L<...<q N mi<-<Mm 



Tr 



Jt(h-B) 



y — y 



AT AT 



f — - ... £ e f« £ detM^ 1 ^) 



(4.51) 



22 Af 

1=0" 1 ^ r,m=0 M=0 nx<—<HM 

x £ [det^SW - i?«) + t) - 1) det N S£{ 

qi<...<q N 

where the tilde over the matrix means, as above, that the temperature factor 
exp{— Y<a=i(ka ~ h + B)/T} is included into the determinant. 

In the thermodynamic limit the numerator of Gj acquire the form 



Tr 



-H/T 



00 I N N-l 

e it{h-B) y< 1 V 
n=o N(N + 1) r m=0p n=0 



2-rri lui 
QN+l Tm N Pn 



xX 



N, P ,r J dq 1 ■ ■ ■ dq N (r] n:m ) (4.52) 

qi<...<q N 

where we denote 

"TV \Vn,m) 

:= detAr^+H^) -H {+ \vn, m )) + (G(x,t) - l)det N S^\rj n>m ). (4.53) 
The function G(x,t) is the thermodynamic limit of the function g(x,t) ( |3.34j ): 



G{x,t) 



2tt 



dke 



-itk +ixk 



(4.54) 



The matrices S^ + '(r] ntm ) and 7o + ) {j} n ,m) have the structure 

8*\*^ = N + 1 + " 2 05 ""'- V, - v 2 , 



(4.55) 
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Matrix elements of the matrices M and Vi,2 (which are N x N dimensional ma- 
trices here) are given by the equations (|4.24|) and Q4.25|) (but with the momenta 
q instead of k). The elements of the matrices 3 (functions TZ\ ^ 3 (q a , q b )) are 

£- h + B E + (q a )E + {q b ) 



^•i +) (?a,gfe) = e" 
K { ^\qa,qb) = e" 



-h+B 1 



47T 



2tt 

£+(g a ) + E + (q b ) 



qi-h+B 



1 



2nE_(q a )E_(q b ) 

The functions _E±, entering ( |4.56| ) are defined in ( f4.26| ), ( ^4.27| ). 

In (|4.52|) , the quantity X N ^ r is the same as in ( |4.20|) above, since 

N 

X N , p , r := £ e^ M £ det M Ufr + \{fi}) 



(4.56) 



M=0 



m<...<^M 



N+l 



6 E det M U^\{\}), 

M=0 \!<...<\ M 



where 



(U^ + \{^})) ab = e^(U^\{^})) ab 

and the elements of the matrices U^ 1,+ ' are given by ( |3.44| ). 
Similarly to ( |4.36| ), we have 



:(+) 

J iV \'ln,m) 



N 



(n 



s *> 



(4.57) 
(4.58) 

(4.59) 



-JV 



where $^ + - ) are some functions which do not depend on n, m. Therefore, the 
result ( |4.39|) for X N)P)r is valid also in the case considered and we can repeat 
the calculation described by equations ( |4.40| )-( |4.42| ), except that we should write 
Hjy instead of S^"^. The last expression in ( |4.42j ) is now written as 



2tt 



dr, 1 + E7" P 



ir)P 



+ e 



- 11 IP 



4 +) (^7)- 



(4.60) 



Here 5^(7757) denotes the expression for S^\r]) given by ft4-53|) in which the 
replacement S^(r]) — > ^yS^^rj) ,1Z^ + \rj) — > r )RS + \rf) has been made. 

Thus we obtain the Fredholm determinant representation for the numerator 
of the correlator 



Tr 



- H / T ^(a;,t)^+(0,0) 



it(h-B)J_ 

2tt 



det ( / + Z + 7 V(v) - 7^ (+) (v) 



+ (G(x,t) - l)det [i + Z + 1 V(r ] ) 



(4.61) 
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where V(r/) is given by ( f4.44| ) and 
1+e cost] ~ (+) 



(+) , 1-g cosr/ ^( +) ^ 



(4.62) 



The integral operators Z, IZ^s possess kernels given by the functions ( |4.15| ), 
(|4.25| ), (|4.56| ) and the function F(j, rj) is given by ( |4.45| ). 

Finally, taking into account the representation fl4.15|) for Tr[e _iT//T ] we obtain 
the following representation for the temperature correlation function 



(4.63) 



where 



G[ + \x,t;h,B) = e u ^- / d V F(j,r]) 

— TV 

det (l + jV(v) - t7^ (+) (^)) + (G(x,t) - 1) det (/ 
V(,) = i±^V 1 + ^V 2 , 

n^( V ) = 1+£ 2 C ° 8V + sinr, HP + l ^ipJL £M (4.64) 

The integral operators Vi j2 possess kernels ( [4.48|) and integral operators 7t^ 3 
possess kernels 

n[ + \q,q') 
n { 2 + \q,q') 

Ki + \q,q') 



E T Aq)E T Aq^ 



1 

47T 



2vr 



El{q>) El(q) 
1 7?(g) ti(q') 



2-k E T {q) El{q'Y (4 ' 65) 

The functions E 1 ^, E 7 are given by equations (^4.49|) with the Fermi weight (|4.50|) . 

Let us now discuss our main results, the representations ( |4.4(j| ) and ( |4.63| ). The 
integral operators V{jf) and IZ^irf) involved into these representations can be 
put in the form usual to the integrable models, i.e., they are "integrable integral 
operators" in the sense of paper (see also ||). To show this we introduce a 
pair of functions 

sin?7 d(k) 



£i( V \k) 



1+e cos 7/ 



Ei(k) + 



2 ~ +v ~ y ' 2 El(k) ' 
The kernels of the operators can be put into the form 
£l(v\k) i T _{k') -F(k) el(v\k') 



t_(k) := Ei(k). (4.66) 



V(7]\k,k') ■■ 

n ( -\k, k') 



n(k-k') 



2tt 



(4.67) 
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This form is important for deriving integrable partial differential equations for 
the correlators and for constructing the corresponding matrix Riemann-Hilbert 
problem. This in turn will make possible the evaluation of different asymptotics 
of the correlators considered. 



5 Particular cases 



Let us discuss now some particular cases of the representations obtained. 

First, consider the "one-component" limit of the theory. Let B — > — oo, 
h — > — oo in such a way that hi = h — B is fixed. Note that h,2 = h + B — > — oo 
in this limit. It means that the energy of the particles of type 2 becomes (see 
(|1.16|) ) infinite and therefore these particles are excluded from the spectrum of the 
theory. So one should has the one-component gas of particles of type 1. Indeed, 
the parameter 7 = 1 in this limit. Remind that the function Fi^f = l,rj) is 
proportional to the 27r-periodic delta- function of the variable rj, and the integrals 
over 77 in (|4.46|) , (|4.63|) can be easily taken. For e = +1 one gets in this way exactly 
the representations for the corresponding temperature correlation functions of the 
impenetrable one-component Bose gas obtained in papers [|L6l O: 



G { f\x,t;h,B) 
G[ + \x,t;h,B) 



-ithi 



h,B — ► — 00 
h-B = h 1 



h,B — > — 00 



det(/ + Vi 



det (J + Vi + K H ) - det (J + V x ) 



(G(x,t) - l)det(/ + Vi 



(5.1) 



where the kernels of the integral operators Vi, TZ^ \ are given by (|4.48| ) and 
( |4.65|) with the Fermi weight 



■&(k) 



1 



k 2 -h 1 

e ~ 



(5.2) 



If e = — 1 (Fermi statistics) then the correlators simplify strongly, reproducing 
the well known simple results for the one-component Fermi gas. The correla- 
tor G% (x,t; h,B) = G[ (x,t; h, —B) vanish in the limit, what means that the 
particles of type 2 are really excluded from the theory. 

Second, consider the particular case of the equal time (t = 0) temperature 
correlators of the two-component gas. It appears that in this case the integrals 
over 77 in ( |4.46 ) and ( 4.63 ) can be also taken with the result 



G[~\x,0;h,B) 



det ( I + 1^-1 v + fi~) 



det ( I + v 



Gi +) (x,0;/i,5) 



det ( / + - — - v + er (+) 



+ (5(x) - 1) det ( I + v 



(5.3) 
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where the integral operators v, f^' possess kernels 

2sin(|x|^^ 



v(k, k') = -yf&{kj 



l -W(jk, k') = ^f&{kj 



n(k - k') 
exp(±ix^- 



2tt 



(5.4) 



with the Fermi weight d{k) given by (|4.50|) . The equal-time correlators (|5.3| ) 
satisfy the relation 



GS +) (x, 0; h, B) = 5(x) + e G\~'(-x 7 0; h, B 



(-)/ 



(5.5) 



which could be also easily derived from the canonical commutation relations 
(|1 . 2| ) . The representation ( |5.3| ) for the equal-time correlation function for the 
impenetrable Fermi gas (e = —1) written in slightly different form has been 
obtained in [TJ[ . 

Discuss now the case of the two-component gas at zero temperature, T = 0. 
In this case one has 



7#(fc) =9(k 2 F -k 2 



(5.6) 



where 9 is the step function and the Fermi momentum is kp = yh-\~ \B\. 

Consider first zero magnetic field B = 0. In this case the ground state is 
degenerate (see Section 1), the parameter 7 = 2, and one obtains 



G[ \x,t;h,0) 



T=0 



-ith 



e "" l — / dri 
2vr ' ' 



G { f\x,t;h,0) 
3 



T=0 



X 



G[ +) {x,t;h,0) 
1 



5 — 4 cos rj 

— /i 

det (J + Vofo) + | n { - ] ) ~ det (/ + V (rj) 
G^ +) (x,t;/i,0) 
3 



T=0 



T=0 



5 — 4 cos rj 



X 



det (/ + V { V ) - n { o"\v)) + t) - 1) det (/ + V ( V ) 



(5.7) 



Here the integral operators act on the functions on the interval (—k F , k F ), e.g., 
(V (t/) ■/)(*) = / dA/Vo(i7|*,A0/(A0- (5-8) 
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The kernels are given as 

£+(v\k) t-{k') - £-(k) £+(vW) 



V ( V \k, k') 

n ( f\k,k') 



n(k - k') 

£-{k)£-{k>) ^ 



2tt 



K +, (v\k,k') 



£ + (r]\k) £ + ( V \k') 



(5.9) 



where 



„ / , , n I + e cos 77 „ . , . sin ri 1 
£ + (rj\k) := '- E + (k) + ' 



7r(l + e cos?]) 
£_(Jfe) := £L(ife), (5.10) 



2 iL(Jfe)' 

and functions E±(k) are given by ( |4.26| ). 

Turn to the case B < 0. In this case the ground state is the Fermi zone filled 
by the particles of type 1; one has 7 = 1 and therefore F(l,r/) = 27rA(r/) where 
A(r/) is a 27r-periodic delta-function. For the correlators one obtains 



_ e -it(h-B) 

G^\x,t;h,B) 

= e it(h-B) 

G { 2 \x,t;h,B) 
G ( 2 + \x,t;h,B) 



B<0,T=0 



det (/ + V (0) + TZ^) - det (/ + D o (0) 



_B<0,T=0 



det (/ + V (0) - 7^ +) (0)) + (G(x,t) - 1) det (/ + V o (0) 
0. 



B<0,T=0 



B<0,T=0 



it{h+B) 



X 



det (J + V (r/) - 7^ +) (r/)) + (G(x, t) - 1) det (/ + V (r/) 



where 



V (0) 



1 + e 



Vi, 



7^ +) (0) 



1 +S S(+) , 1 ~ e i3(+) 



-^1 



/c 3,0 • 



(5.11) 



(5.12) 



The integral operators act on the functions on the interval (—kp, kp) by the rule 
(|5.8|). The kernels are given by 



Vi, (MO 

^-i,o(k, k ) 



E + (k)E_(k') - E_(k)E+(k') 

n(k - V) 
E + (k)E + (k') 



2tt 



^•3,0 {k, k ) 



2nE^(k)E_(k')' 



(5.13) 



and by (p|). 

If B > 0, the ground state is formed by the particles of type 2, and the 
formulae for the correlators are obtained in an obvious way, due to the relation 



34 



Conclusion 



In this paper the determinant representation for the two-point time-dependent 
temperature correlation functions of the one-dimensional two-component bosons 
and fermions are obtained. The Fredholm determinant of the integrable integral 
operator enters our answers ( 4.46|) , ( |4.63|) . It makes possible deriving intergrable 



differential equations for the correlators and calculating their large time and 
distance asymptotics. 

From the technical point of view, the essential thing in our approach is using 
the XXO eigenfunctions for the two-component model with the infinitely strong 
coupling. This gives the opportunity of explicit calculation of the correlators. The 
same technique can be applied for calculating temperature correlation functions 
in the Hubbard model in the infinitely strong coupling limit. We are going to 
present the corresponding results in the next publication (with N. I. Abarenkova). 
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Appendix 

Here we consider matrix elements of the matrices , , entering the 
representations for normalized mean values for the model on the finite interval 
(Section 3). Our aim is to rewrite matrix elements of these matrices in terms of 
functions which are well defined in the thermodynamic limit. We prove also the 
representations (|3.12|) , (|3.36|) for the matrices considered. 

We begin from a pair of indentities 

cot7TZ = - + £ .__ I t ^- (A.l) 

z \z - irj 7tj J sin tiz j^ooiz-nj) 2 

The value of the r.h.s. do not depend on organizing of summation. Let us fix the 
method of summation in such a way that the term 1/irj does not contribute to 
the sum in the first equation, e.g., combining terms with the same \j\ together 
and summing up in increasing order of \j\ = 0, 1, 2, . . .. The sums below should 
be understood with this prescription (where necessary). 
In terms of quasimomenta k and q, 

. . 2tt / A+e^N - 1 A . _ . . n . 

(*)i = T (-(— )— + ') + L' JGZ ' (A ' 2) 
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the both indentities ( |A.l ) become (g = q + A/L): 

11 + euju\ 



2 ^ 1 

Ly q-k 



tan 



V 2 



L 2 



E 



r' (g - k)< 



(A.3) 



where u; := expzA and v := expiB. The expression below, due to the first 
indentity in (|A.3|), can be rewritten as 



11 + eijjv 



,-,,,12 



E 



-itq 2 +ixq 



11 + euv\ 



^—itq 2 +ixq ^—itk 2 +ixk ^—itk 2 +ixk 



q — k 



+ E 



q — k 



= i + ^os(A-e) eH(fc) + eS in(A-e) [e _ (A;)] _ 2 

where the functions e_(fc) and e^~^(fc) are introduced, 



„(-) 



g—itq 2 +ixq ^—itk 2 +ixk 



q — k 



6 



e_(/c) = exp 



'itfc 2 — ixk 



(A.4) 



(A.5) 



Let us consider the elements of the matrix given by 

„ — itq 2 +ixq 



(S^) ab = e4k a ) e^h) 



|1 + EVUj\ 2 

T 2 



E 



(q - k a )(q - h) 



(A.6) 



Consider first the off-diagonal elements, i.e., the case when a ^ b, hence k a ^ k b . 
For the off-diagonal elements of matrix S^~\ due to ( |A.4|) , one has 



(S^) ab = e.(k a ) e.(k b ) 



L 2 (k a - k b 



■E 



^—itq 2 +ixq ^—itq 2 +ixq 



q-k a q-h 



2 / l + e cos(A - 9) e ( + \k a )e_{k b ) - e_(fc a )ej \k b ) 



e sin(A - 9) [e_(A; a )]^ 1 e-(fc b ) - e-(k a )[e-(k b )]^ 1 
2 k a k b 



(A.7) 



Let us consider now the diagonal elements of the matrix '. Using first the 
second indentity in (|A.3|) and then ( |A.4j ), one has for the diagonal elements 



I n n C- 



|1 + euv\ 



■ E 



e -itq 2 +ixq _ g 



{q - k a ] 



+ e 
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2 f 1 + e cos(A - 0) 
~ L\ 2 

2 e"^ 2 ^ - e -^+^(l + i(x - 2tfc a )(g - k a 

L 2 ? (q - k a ) 2 

e sin(A - 6) . , . ] 
+ i(x - 2tk a ) + 1. 



(A.8) 



It is easy to see that the expressions in (|A.8| ) following 1 + e cos(A — 6) and 
esin(A — B) are the limiting cases of the corresponding expressions in ( |A.7| ) when 
K — > k a . Thus, collecting ( |A.7|) together with flA-ST) , one can conclude that the 
matrix has the structure (|3.12|) where matrices V}^ are given by (|3.13| ) with 
the diagonal elements being understood in the sense of the l'Hopitale rule. 

Consider now matrix elements of the matrices and R( + > given by ( |3.35|) . 
In this case, instead of ( |A.3| ), the following indentities are useful 



-V — 

L^k-q 

k 



—e 



M-0 

tan 

V 2 



|1 + eujv\ 

T 2 



E 



(k-qf 



(A.9) 



which are, of course, direct consequences of (|A.1|) . Also, instead of ( |A.4|) , one has 



EU)h>\ 2 



L 



E 



—itk 2 +ixk 

k — q 



l + e cos(A - 6) 



,(+) 



e sin(A - 6) 



-(?)]• 



(A.10) 



where the function e^ + \q) is defined as 

Vixk p —itq 2 +ixq 



,(+) 



L ^ 

k 



-itk 2 



k — q 



(A.n; 



For the matrix S^ + ' all the calculation explained above on the example of the 
matrix can be done similarly This leads to the structure ( [3.36| ) of the 
matrix S^ + '. The matrix defined by ( |3.35|) can be put into the form (|3.36| ) 



by means of the indentity (|A. 10|) applied to each factor defining R ab in ( |3.35| ) . 

Let us now consider the thermodynamic properties of the functions containing 
sums over quasimomenta q (or k) in their definitions. These functions are 
entering the matrix elements of the matrices R^ and the function g(x,t) 
entering the representations (|3.45| ), ( |3.46| ). 

Let us begin from the simplest example, the function g, which is defined as 



g(x,t) 



k 



-itk 2 +ixk 



k = k + 



A 



(A.12) 
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In the limit L — > oo the values (k)j of the quasimomenta k fill the interval 
(—00, +00) densely, (k)j + i — (k)j = (k)j + i — (k)j = 2n/L, and the sum over 
values of k should be changed for the integral, 



00 - 00 



x£/(*o : = 1 e /((%) -> ^ y <**/(*)• (a. 13) 

fc J=-oo -OO 

Thus, in the limit L — > 00 for the function g one has 



c/(x,t) -> G(z,t) = — / dke~ itk ' 2+ixk . (A.14) 

—00 

The dependence on A of the function g disappears in the limit L — > 00. The 
same happens with respect to the dependence on m of functions g m := g| A= 2gm . 
All functions <? m have the same thermodynamic limit which is equal to function 
G. 

Consider the function eW(g) given by ( |A.11| ). Recall that the sum is under- 



stood in the sense of prescription described above. This prescription regularizes 
the sum on the "infinities" while the subtraction term, in the limit L — » 00, re- 
move the singularity coming from the pole on the integration counter. It means 
that such a sum turns into the integral in the sense of the principal value, i.e., 
for the function e' + '(g) one has 

e i - + \q)^E{q) = P.v. dk— -. (A.15) 

J 7i(k — q) 

—00 

Note that the dependence on A and of the function e^ + '(q) disappears in 
the limit L — > 00 (function E(q) does not depend on A and O) since k and q 
(function's argument) are assumed to run the interval (—00, +00) continuously. 

The function e^~\k) given by ( |A.5|) has the same thermodynamic limit as the 
function e <y+ \k): 



e 



00 

(-) 



dq— — . (A.16) 

n(q — k) 



Thus, all functions determining matrix elements of the matrices and 
are well defined in the thermodynamic limit. 
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